Generalized Coarse Graining Procedures for Flow in Porous Media     

Sabine Attinger 《Computational Geosciences》2003,7(4):253-273

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows     

Frederico Furtado  Felipe Pereira 《Computational Geosciences》2003,7(2):115-135

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.  相似文献   

地震波研究://car          下载免费PDF全文

张海明  陈晓非 《地震学报》2003,25(5):465-474

回顾了近4年来中国地球物理学家在地震波研究方面取得的主要进展. 分别简要地讨论了层状均匀介质、横向不均匀介质和多孔介质中的地震波传播，地震面波及波形反演，以及地震勘探、测井问题中的地震波研究方面取得的新成果. 指出当前地震波研究的重要课题是发展高效的数值方法，并利用该方法开展复杂介质中的地震波激发和传播，以及强地面运动的研究，为精细的地震危险分析与预测奠定基础.   相似文献   

Confined aquifer as wave-guide and its responses to geo-acoustic waves     

简文彬  陈葆仁  卢华复 《地震学报(英文版)》1997,10(3):387-395

ＣｏｎｆｉｎｅｄａｑｕｉｆｅｒａｓｗａｖｅｇｕｉｄｅａｎｄｉｔｓｒｅｓｐｏｎｓｅｓｔｏｇｅｏａｃｏｕｓｔｉｃｗａｖｅｓＷＥＮＢＩＮＪＩＡＮ１）（简文彬）ＢＡＯＲＥＮＣＨＥＮ２）（陈葆仁）ＨＵＡＦＵＬＵ２）（卢华复）１）Ｄｅｐａｒｔｍｅｎｔ...  相似文献   

岩体水力学基础(三)──岩体渗流场与应力场耦合的集中参数模型及连续介质模型     

仵彦卿 《水文地质工程地质》1997,(2)

本文在分析了岩体系统的结构性和水力学特征后，提出了岩体渗流场与应力场耦合数学模型的机理分析、混合分析及系统辨识建模方法。运用系统辨识方法建立了岩体渗流场与应力场耦合的集中参数模型，并应用于解决实际问题；运用机理分析和混合分析方法建立了岩体渗流场与应力场双场耦合及与温度场三场耦合的连续介质分布参数模型。  相似文献   

Inferring the relative measure of principal stress components     

Peiliang Xu  Seiichi Shimada 《Geophysical Journal International》1997,130(1):F1-F4

The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.  相似文献   

Wide-angle seismic velocities in heterogeneous crust   总被引：5，自引：0，他引：5  

John Brittan  Mike Warner 《Geophysical Journal International》1997,129(2):269-280

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.  相似文献   

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation   总被引：3，自引：0，他引：3  

Vadim Levin  Jeffrey Park 《Geophysical Journal International》1997,131(2):253-266

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.  相似文献   

渗流场反分析及其在多孔介质中的应用     

许模 《地质灾害与环境保护》1996,(2)

本文讨论了复杂地质环境中渗流场的反分析问题，利用地下水流数值模拟结合最优化技术，通过对地下渗流场的模拟和寻优，求得描述地下含水介质非稳定渗流场的重要参位──渗透系数和给水度的模型解，并应用于各向同性的多孔介质渗流场反分析实例中。  相似文献   

多媒体创作工具Compel的应用     

徐晓平  陈天与 《物探化探计算技术》1996,18(4):325-330

本文介绍了一种多媒体应用的创作工具。利用这种工具可以编辑和组织文字、图象、图表、声音以及视频等多种媒体元素。并通过媒体链接以及超级链接，完成多媒体应用的制作。同时，也介绍了多媒体应用的播放。  相似文献